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Abstract. The effect of environmental temperature on neuronal spiking behaviours 
is investigated by numerically simulating the temperature dependence of spiking 
threshold of the Hodgkin-Huxley neuron subject to synaptic stimulus. We find that 
the spiking threshold exhibits a global minimum in a "comfortable temperature" range 
where spike initiation needs weakest synaptic strength, indicating the occurrence 
of optimal use of synaptic transmission in neural system. We further explore the 
biophysical origin of this phenomenon in ion channel gating kinetics and also discuss 
its possible biological relevance in information processing in neural systems. 
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Neuronal spiking is the basis of communication and coding in nervous systems 
and has been studied for decades [11 [2l |3]. Spiking behaviours of a neuron come from 
the activation and inactivation of various ion channels, which may vary due to the 
influences of environmental factors [H [5]. Temperature is one of the most important 
regulators to neuronal activities [SI [71 El E] . Electrophysiological experiments on neural 
spiking activities are mostly conducted either at room temperature in vitro or at body 
temperature in vivo, thus our knowledge about how a neuron behaves in a continuously 
changing temperature environment is usually fragmental. Variation in temperature 
leads to two main impacts in the spiking dynamics of excitable neurons: one is on 
the maximum ion channel conductances [TDl E] and the other is on ion channel gating 
kinetics [121 [ISl E] ) hence altering the shape and amplitude of action potentials [15] and 
the generation [16l [17] and propagation [18] of spikes. Neuronal spiking activities can 
be described by means of firing properties, such as firing rates [l9j, the precise timing of 
spikes [201 [211 [22], or the spiking threshold behaviors [231 1211 [251 126]. The temperature 
dependence of spiking threshold of squid giant axons under current injection has been 
investigated experimentally [221 [211 I2Z] and numerically [211 |29]. More realistically, 
in fact, neurons communicate with each other via synaptic connections and respond 
to synaptic stimuli by firing spikes, so as to detect, transmit, and process neural 
information. It is an interesting question how the environmental temperature affects 
the spiking properties of a neuron stimulated by realistic synaptic input. 

This study is to approach the above topic by examining the temperature dependence 
of spiking threshold of a neuron subject to synaptic stimulus, based on widely-accepted 
Hodgkin-Huxley (HH) neuron model, as an example, which was originally proposed 
to account for the excitable properties of giant squid axons [12] and is now regarded 
to be an useful paradigm that accounts naturally for the spiking behaviors of real 
neurons. Our results reveal that there is an environmental temperature range where 
the spiking threshold shows a global minimum, implying that the optimal use of 
synaptic transmission occurs in neural systems. The biophysical origin of this interesting 
phenomenon together with its biological relevance is also investigated and discussed. 

The dynamics of the HH model with synaptic stimulus is described by the following 
coupled differential equations: 

dV/dt = {I syn - I ion ) I C, 

dmjdt = [m^{V) -m]/TrniV,T), 

dh/dt = [h^iV)-h]/ThiV,T), ^ ^ 

dn/dt = [n^{y)-n]/Tn{y,T), 
with V being the membrane potential, C the membrane capacity, n the activation 
variable of potassium channel, and m and h the activation and inactivation variables 
of sodium channel, respectively, moo, ^oo, ""-oo and r^, t^,, t„ represent the saturated 
values and the time constants of the gating variables, respectively. The r's in Eqs. ([1]) 
depend on the environmental temperature T, usually by being divided by a Qio factor 
of a, Qio{T, a) = 0*^"^"^°)/^° with the a value being usually chosen as a = 3 (suggested in 
Ref. [1]) and Tq = 6.3 °C denoting the experimental temperature for the original model 
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construction [12] . The ionic current lion includes the usual sodium, potassium, and 
leak currents: 

hon = GNa{T)m^h{V - VNa) + GK{T)n\V - Vk) + Gl{T){V - Vl)^ (2) 

where V/va, Vl are the reversal potentials for the channel currents. The maximum 
channel conductances Giya, Gk, and Gl are regulated by environmental temperature 
with a Qio factor of a = 1 — 1.5 (see Ref. [1]). The synaptic input is modeled by 
hyn = gsynit)(y — Vgyn) with Vgyn being the synaptic reversal potential and Qsynit) 
being the time-dependent post-synaptic conductance, Qsynif) = GsynOiit — to), where 
Iq represents the onset time of the synapse, Ggyn determines the peak of synaptic 
conductance and a{t) = {t / Tsyn) exp{l — t/Tsyn), t > 0, with Tgyn determining the 
characteristic time of the synaptic interaction. In this study we choose Tgyn = 2 ms, 
and Vsyn = mV to mimic excitatory synapse input in neural system [251 EO] • The 
other values of parameters can be found in Ref. [12]. The spiking threshold here is 
characterized as the critical value of Ggyn, by which the membrane potential of the 
stimulated neuron exceeds a voltage threshold Vth (chosen as Vth = — 20mV here). 

Our investigation begins with the question how environmental temperature affects 
the spiking threshold of HH neuron stimulated by synaptic input. In figured! one may 
find that the spiking threshold undergoes a U-shaped dependence on the environmental 
temperature, i.e., there is a global minimum spiking threshold in a temperature range. 
We refer to this temperature range as "comfortable temperature" one for the neuron. In 
the comfortable temperature range the neuron needs weakest synaptic stimulus to fire 
spikes; from the engineering perspective this facilitates the extraction, transmission and 
processing of information in neural systems. Our result suggests that environmental 
temperature can make the neurons communicate easier with each other by maximizing 
the utility of synaptic transmitters. The optimal use of synaptic stimulus in nervous 
system is of great biological significance for real neurons. 

Theoretically, a U-shaped dependence of a physical quantity characteristic may 
result from the competition of (at least) two contrary factors/aspects. In a neural 
system, environmental temperature influences the maximum ion conductances and the 
channel gating rates, as demonstrated by experimental observations [4j and explicitly 
included in our model simulation. To find out what factors account for this U- 
shaped dependence phenomenon, we make a comparative examination. The effects 
of temperature regulation on spiking threshold through the maximum ion conductances 
and through the gating time constants of channel variables (see figure [2]) are investigated 
respectively. If only the influence of temperature on the maximum ion conductances is 
considered, the spiking threshold increases monotonously with temperature. In direct 
contrast, the sole temperature effect via gating kinetics of ion channels gives rise to a 
temperature dependence of spiking threshold that almost resemble the control result in 
figure [H Thus one can reach the conclusion that it is the dominant role of temperature on 
the gating kinetics of ion channels that yields the phenomenon of optimal use of synapse 
transmission in neural system, though changes in maximum channel conductances 
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simultaneously contribute slightly to the elevation of spiking threshold with increasing 
environmental temperature. 

So far we have been aware of that the occurrence of optimal use of synaptic 
transmission attributes mainly to the thermal impacts on the gating kinetics of ion 
channels, the direct link lacks between these two: why does a monotonous dependence 
of gating time constants on temperature result in a non-monotonous temperature 
dependence of spiking threshold in neural systems? As well known, the activation of 
sodium ion channel depolarizes the membrane potential and forms the rise phase of an 
action potential; on the contrary, the activation of potassium ion channel along with 
the inactivation of sodium ion channel forms the decay phase of an action potential. 
Is it the competition of the gating behaviours between these two types of channels 
that yields the U-shaped temperature dependence of spiking threshold? The answer to 
this question is unfortunately negative (simulation results are not shown here). Thus 
we divide the three gating variables of ion channels into two categories: one is the 
sodium ion channel activation variable m, which is beneficial for a neuron to initiate 
spikes, and the other includes the sodium ion channel inactivation variable h and the 
potassium ion channel activation variable n, which serve to terminate action potentials. 
We speculate that it is the interaction between these two contrary (competitive) aspects 
that lead to the temperature dependence of the neuron's spiking threshold. To confirm 
our interpretation, we design further comparative explorations. In figure [3] one can find 
contrasting dependencies of spiking threshold on environmental temperature between 
the thermal impact via activation kinetic of sodium ion channel alone and that via the 
inactivation of sodium ion channel together with the activation kinetic of potassium ion 
channel. These results are in good agreement with our speculation. 

In the above simulation results for the temperature dependence of spiking threshold 
are obtained by choosing the specific Qiq factors of a = 3 for time constants of channel 
gating variables and a = 1.25 for the maximum ion conductances. Although other Qio 
factors might also be used for the neural system [U [31] , they are not expected to yield 
essential difference and the U-shaped temperature dependence still occurs. Examples for 
the HH neuron are shown in figure HI where one can see that temperature dependencies 
of spiking threshold are very similar to the control case (see figure [1]). The characteristic 
time Tsyn of synaptic interaction may also vary with environmental temperature, but 
again it does not lead to any essentially different result against our conclusions. In 
fact, the elegance of optimal use of synaptic transmission in a range of environmental 
temperature in squid giant axons can be further generalized to other nervous systems, 
such as the cochlear nucleus neuron in auditory system [32]. Presumably, this subcellular 
mechanism of thermal impacts on the neuronal spiking via ion channel kinetic behaviours 
can serve as the base for the psychophysical evidence that various animals live optimally 
in a specific range of environmental temperature. 

In summary, we have investigated the effects of environmental temperature on the 
spiking behaviours of a neuron subject to synaptic stimulus. Our results reveal that 
the spiking threshold shows a global minimum in the so-called comfortable temperature 
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range, implying the occurrence of optimal use of synaptic transmission in neural system. 
We further illustrate that the emergence of this phenomenon attributes mainly to the 
combined competition of the temperature-dependent gating kinetics of ion channel 
activation/inactivation variables. The optimal use of synaptic transmission in neural 
system can, from an engineer's perspective, largely facilitate the extraction, processing 
and transmission of neural information in neurons. 
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Figure 1. The dependence of spiking threshold of an HH neuron subject to synaptic 
stimulus on environmental temperature. Qio factors of a = 1.25 and a — i are 
employed for the maximum ion conductances (the G's in Eq. ([2|)) and the gating time 
constants of ion channels (the t's in Eqs. ([T])), respectively. 




Figure 2. The influence of temperature on the spiking threshold via maximum channel 
conductances and via channel gating kinetics , respectively. Dots: only the maximum 
ion channel conductances (the G's in Eq. ([2])) are multiplied by a Qiq factor of a = 1.25. 
Stars: only the time constants of gating variables (the t's in Eqs. ^) are divided by 
a Qio factor of a = 3. 
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Figure 3. The influence of temperature on tlie spiicing thresiiold via the thermal 
regulations on the activation/inactivation variables. If only the time constant of sodium 
channel activation variable m is regulated with a reciprocal Qio factor of a = 3, the 
spiking threshold decreases with increasing temperature (triangles); if only the time 
constants of sodium channel inactivation variable h and potassium channel activation 
variable n are divided by a Qio factor of a = 3, spiking threshold increases with 
temperature (dots). 
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Figure 4. The temperature dependencies of spiking threshold of the HH neuron in 
a variety of Qio factors for maximum ion conductances (the G's in Eq. ([2])) and the 
gating time constants of channel variables (the t's in Eqs. ([1])). The Qiq factors are 
chosen within the physiologically suitable range. 



